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1. — Introduction 



For many years it has been known that locahzed (non-dispersive) solutions exist 
to the (homogeneous) wave equation[l,2,3], endowed with subluminal or Superlumi- 
nal[4,5,6,7,8,9] velocities. These solutions propagate without distortion for long distances 
in vacuum. 

Particular attention has been paid to the localized Superluminal solutions like the 
so-called X- waves [5, 6, 8] and their finite energy generahzations[7,8]. It is well known that 
such Superluminal Localized Solutions (SLS) have been experimentally produced in acous- 
tics[10], optics[ll] and more recently microwave physics[12]. 

As is well known, the standard X-wave has a broad band frequency spectrum, starting 
from zero[8,9] (it being therefore appropriate for low frequency applications). This fact 
can be viewed as a problem, because it is difficult or even impossible to define a carrier 
frequency for that solution, as well as to use it in high frequency applications. 

Therefore, it would be very interesting to obtain exact SLSs to the wave equations 
with spectra localized at higher (arbitrary) frequencies and with an adjustable bandwidth 
(in other words, with a well defined carrier frequency). 

To the best of our knowledge, only two attempts were made in this direction: one by 
Zamboni-Rached et a/. [8], and the other by Saari[13]. The former showed how to shift 
the spectrum to higher frequencies without dealing with its bandwidth, while the latter 
worked out an analytical approximation to optical pulses only. 

In this work we are presenting analytical and exact Superluminal localized solutions 
in vacuum, whose spectra can be localized inside any range of frequency with adjustable 
bandwidths, and therefore with the possibility of choosing a well defined carrier frequency. 
In this way, we can get (without any approximation) radio, microwave, optical, etc., 
localized Superluminal waves. 

Taking advantage of our methodology, we obtain the first analytical approximations 
to the SLS's in dispersive media (i.e., in media with a frequency dependent refractive 
index) . 

One of the interesting points of this work, let us stress, is that all results are obtained 
from simple mathematical operations on the standard "X-wave" . 
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2. — Super luminal localized waves in dispersionless media 

Let us start by dealing with SLSs in dispersionless media. From the axially symmetric 
solution to the wave equation in vacuum (n = 1), in cylindrical coordinates, one can easily 
find that 

ij{p,z,t) = Jo(Me+''^'e-^"* (1) 

with the conditions 

^2 = ^-^.^ fc2>o, (2) 

where Jo is the zeroth-order ordinary Bessel function; and kp are the axial and the 
transverse wavenumber respectively, uj is the angular frequency and c is the light velocity. 

It is essential to call attention right now to the dispersion relation (H). Positive 
(but not constant, a priori) values of A;^, with real kz-, do allow both subluminal and 
Superluminal solutions, while implying truly propagating waves only (with exclusion of 
the evanescent ones). We shall pay attention in this paper to the Superluminal solutions. 
Conditions (^ correspond in the {uj,kz) plane to confining ourselves to the sector shown 
in Fig.l; that is, to the region delimited by the straight lines uj = ±ckz- 

Important consequences can be inferred from (0), when performing the coordinate 
transformation 

I = ^ (3) 

I kz = (uj/c) cos 6 , 

which yields 

kp = (lu/c) sin 6 . (4) 
With this transformation, solution (|l|) can be rewritten, in the new coordinates {uj,9), 

as 

HpX) = Jo(-psin^)e+^-^™^^ , (5) 
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where C, = z — Vt, and where the propagation speed (group velocity) is obviously V = 
cj cos Q. Equation (|^) states that the beam is transversally localized in energy, and 
propagates without suffering any dispersion. It should be noticed also the relationship 
between V and Q: namely, each value of Q yields a different wave velocity. This fact will 
be used in the next Section. 

Equation (^) represents the well known "Bessel beam". As can be seen, such an 
equation has two free parameters, uj and Q. Considering Q constant, and making a su- 
perposition of waves for different frequencies, one can obtain localized (non-dispersive), 
Superluminal solutions; namely 

^(p,() = / 5(cj)Jo(-psine) e+'"^™^Mcu . (6) 
h c 

In eq.(^, if an exponential spectrumB like S{uj) = e~"'^ is considered, one obtains, by 
use of identity (6.611.1) of ref.[14], the ordinary X-shaped wave: 



where a is a positive constant. 

This solution is a wave that propagates in free space without distortion and with the 
Superluminal velocity V = 1/ cos 0. Because of its non-dispersive properties, and its low 
frequency spectrum, the X-wave is being particularly applied in fields like acoustics[5]. 
The illustration of an X-wave, with parameters a = 10~^ s and = 5 c, is shown in Fig.2. 



3. — Superluminal localized waves for arbitrary frequencies and 
adjustable bandwidths 

In the last Section, it has been shown that a superposition of Bessel beams can be 
used to obtain a localized and Superluminal solution to the wave equation in a disper- 
sionless medium. It is known that it may be a difficult task, it being possible, or not, 
finding analytical expressions for eq.(^). Its numerical solutions usually brings in some 

*It is easy to see that this spectrum starts from zero, it being suitable for low frequency applications, 
and has the bandwidth Ato — 1/a 
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inconveniences for further analysis, uncertainties concerning the fast oscillating field com- 
ponents, etc.; besides implying a loss in the physical interpretation of the results. Thus, 
it is always worth looking for analytical expressions. 

Actually, the kind of solution found by us for eq.(|^) is strictly related to the chosen 
spectrum S{uj). Following previous work of ours[8], we are going to present our spectrum 
together with its main characteristics. 

3.1 — The S{u) spectrum 

One of our main objectives is finding out a spectrum which can preserve the integra- 
bility of eq.(|^) for any frequency range. In order to be able to shift our spectrum towards 
the desired frequency, let us locate it around a central frequency, with an arbitrary 
bandwidth Au. 

Then, let us choose the spectrum 

5H = (^)'"e— (8) 

where V is the wave velocity, while m and a are free parameters. For m = 0, it is 
S{lj) = exp[—auj], and one gets the (standard) X-wave spectrum. 

After some mathematical manipulations, one can easily find the following relations, 
valid for m ^ 0: 

= — ; ^ ; ; — ; rr (8.1) 

(Acu±/cuc) - In (1 + (Acui/wc)) 

TTl 

00, = -. (8.2) 
a 

Here, because of the non-symmetric character of spectrum (^), let us call Adj^ (> 0) 
the bandwidth to the right, and Au;_ (< 0) the bandwidth to the left of ujc\ so that 
Acj = Au;+ — Acj-. It should be noted however that, already for small values of m 
(typically, for m > 10), one has Ac<j+ —Auj-. Once defined Uc and Au, one can 
determine m from the first equation. Then, using the second one, a is found. 

Figure 3 illustrates the behavior of relation (8.1). From this figure, one can observe 
that the smaller Auo/uOf. is, the higher m must be. Thus, one can notice that m plays the 
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fundamental role of controlling the spectrum bandwidth. 

From the X-wave spectrum, it is known that a is related to the (negative) slope of 
the spectrum. Contrarily to a, quantity m has the effect of rising the spectrum. In this 
way, one parameter compensates for the other, producing the localization of the spectrum 
inside a certain frequency range. At the same time, this fact also explains (because of 
relation (8.2)) why an increase of both m and a is necessary to keep the same Uc- This 
can be seen from Fig. 4. 

In Fig. 4, both spectra have the same Uc. Taking the narrow spectrum as a reference, 
one can observe that, to get such a result, both quantities m and a have to increase. 
Moreover, this figure shows the important role of m for generating a wider, or narrower, 
spectrum. 

3.2 — X-type waves in a dispersionless medium 

To illustrate the use of the proposed solutions, let us define the ordinary X-wave, by 
rewriting eq.(|^) with S{ijj) = exp[— atu]: 



^(P,C) = /o°°^o(^PV'^g^-l)e-('^^-<)-/^d.; 

= ^/o°°^o(^Pv'^^§?^)e-("^-<)-/M^) 

which is the same as 



(9) 



X = ^(p,C) = ^ = , (10) 

where Uq is the refractive index of the medium underlying these considerations. Applying 
our spectrum expressed by eq.(^, equation (^ can be rewritten as 



*(P.C) = V j^('^rj„ {^p{n^ - 1 e-<-'- 0-/' d(^) , (11) 
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We have therefore seen that the use of a spectrum hke (^) allows shifting it towards 
any frequency and confining it within the desired frequency range. In fact, this is one of 
its most important characteristics. 

It can be seen that eqs.(^ and (^) are equivalent. Deriving eq.(|^) with respect 
to {aV — iC), a multiplicative factor (y) is each time produced (an obvious property of 
Laplace transforms). In this way, it is possible to write eq.(|lT]) as : 

A different expression for eg . ([TT|) , without any need of calculating the m-th derivative 
of the X-wave, can be found by using identity (6.621) of ref. [14]: 



T/ ^^ r(m + l)X"*+i r./rn+l m ^ , ^ yv \ 




where X is the ordinary X-wave given in eq.(0), and F is a Gauss' hypergeometric 
function. Equation (12') can be useful in the cases of large values of m. 

Let us call attention to equations ([T2|) and (12'): to our knowledge,! no analytical 



expression had been previously met for X-type waves, which can be localized in the neigh- 
bourhood of any chosen frequency with an adjustable bandwidth. Equations (0) allow 
getting one or more of them in a simple way: All that has to be done is calculating the 
m-th derivative of X with respect to [aV — Alternatively, one can have recourse to 
eq.(12'). 

Fig. 5 shows an example of an X-shaped wave for microwave frequencies. To that 
aim, it was chosen tUc = 6 x 10^ GHz and Auj = 0.9 oOc, and the values of n and a were 
calculated by using eqs.(8.1) and (8.2): thus obtaining m = 10 and a = 1.6667 x 10^^. 
As one can see, the resulting wave has really the same shape and the same properties as 
the classical X-waves: namely, both a longitudinal and a transverse localization. 



'I'lt can be noticed that dX/d{aV ~ iC) — {iV)^^dX/dt. Time derivatives of the X-wave have been 
actually considered by J.Fagerholm ct al.[15]: however the properties of the spectrum generating those 
solutions (like its shifting in frequency and its bandwidth) did not find room in that previous work. 
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4. — Superluminal localized waves in dispersive media 



We shall now pass to dealing with dispersive media. 

In Section (2), equations (|]) and (H) were written for a dispersionless medium 
[riQ = constant, independent of the frequency). However, for a typical medium, when 
the refractive index depends on the wave frequency, nlu), those equations become[13] 

(13) 

k^{u) = ^n{uo) cos{9) . 

The above equations describe one of the basic points of this work. In Section (2) it 
was mentioned that 6 determines the wave velocity: a fact that can be exploited when one 
looks for a localized wave that does not suffer dispersion. In other words, one can choose 
a particular frequency dependence of 6 to compensate for the (geometrical) dispersion 
due to the variation with the frequency of the refractive index[13]. 

If the frequency dependence of the refractive index in a medium is known, within a 
certain frequency range, let us see how the consequent dispersion can be compensated 
for. When a dispersionless pulse is desired, the constraint = a + u b must be satisfied. 
And, by using the last term in eq. (p^ , one infers that such a constraint is forwarded by 
the following relationship between 6 and u: 

cos{9iu)) = , 14 

uj n[uj) 

where a and b are arbitrary constants (and b is related to the wave velocity: b = 1/V). 
For convenience, we shall consider a = 0. Then, eq.(B|) can be rewritten as 



vI/(p,C) = ^"5(c.)Jo(pc.^^-62)e+-^^d^. (15) 

Let us stress that this equation is a priori suited for many kinds of applications. In fact, 
whatever its frequency be (in the optical, acoustic, microwave,... range), it constitutes the 
integral formula representing a wave which propagates without dispersion in a dispersive 
medium. 



8 



Now, let us mention how it is possible to realize relation (|Tip for optical frequencies. 
Although limited to the case of the air, or of low- dispersion media, the axicon[5, 6, 11, 16] is 
one of the simplest means to realize it. Another possibility is using "spectral hole burning 
filters", or holograms [17]. More in general, one can follow a procedure similar to the one 
illustrated in Figure 6. 

The process illustrated in Fig. 6 is actually simple. In fact, there is a different deviation 
of the wave vector for each spectral component in passing through the chosen device 
(axicon, hologram, and so on): and such a deviation, associated with the dispersion due 
to the medium, makes the phase velocity equal for each frequency. This corresponds to 
no dispersion for the group-velocity. More details about the physics under consideration 
can be found in Ref.[13]. 

Now, let us consider a nearly gaussian spectrum as that given by eq.(|), and assume 
the presence of a dispersive medium whose refractive index (for the frequency range of 
interest) can be written in the form 



n{uj) = no + uj 6 , (16) 

where no is a constant, while 5 is a free parameter that makes it possible a linear behavior 
of n{uj): something that is actually realizable for frequencies far from the resonances asso- 
ciated with the used material. Notice that the linear relationship between the refractive 
index and the wave frequency assumed in eq.(|16|) is not necessary: but its existence gets 
our calculations simplified. 

In this way, substituting eq.([T^) into eq.(|I^), and considering the spectrum, shifted 



towards optical frequencies, given by eq.(||), a relation similar to eq.(|T^ is found 



npX,S) = V /f (^)Vo (^^P\/^K + M^-lj e-(^^-^WM^) • (17) 
To the purpose of evaluating eq.(0), let us make a Taylor expansion and rewrite it 



as 



vl>(p,C,5) = ^{p,C,0) + ^ ^\s=o + -Q^\s=o + ^ -^\s=o + ... (18) 

For the above equation it is known that, if 6 is small enough, it is possible to truncate 
the series at its first derivative. For the time being, let us assume this is the case and 
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that there is no problem on truncating eg. (P^ . One can check Fig. 7, which shows typical 
values of 6 for Si02, a typical raw- material in fiber optics. 

Looking at eq.(|l^), one can notice that its first term \E'(p, ^, 0) is already known to us, 
because it coincides with the solution given by our eg. (p!T]) . To complete the expansion 
(0), one must find ^\s=q. After some simple mathematical manipulations, one gets 
that 



93 c2 Jnl^^ - 



(19) 

This integral can be easily evaluated by using identity 6.621-4 of ref.[14], so to obtain 



W l« - (-1'"'" awzno^. - '0 ^1 . (20) 

As in the case of eqs.(12), (12'), another ioTm for expressing eq.(|19|) can be found by 
having recourse once more to the identity (6.621) of ref.[14]: 



(9^ _ nop2r(m + 4)X"^+4 /m + 4 -m - 1 n 2 



where, as before, X is the ordinary X-wave given by eg. ([10|) and F is again a Gauss' 
hypergeometric function. Once more, equation (20') can be useful in the cases of large 
values of m. 

Finally, from our basic solution ([T2|) and its first derivative (pOD , one can write the 
desired solution of eq. ([T5|) as 



an V T/3 „ fin+2 

(21) 
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However, if one wants to use equations (12') and (20'), instead of eqs.(|T2p and (PP]), 
the solution ( pi]) can be written in the form 



r(m + l)X"^+^ j-.(m+l m ^V'^ ^ 2^^^ 



(21') 



nop2r(m+4)X"^+4 ( m+i -m-l . q . ^^2 i ^ ^2 A 



It is also interesting to notice that, e.g., the approximated Superluminal localized solu- 
tion (pID for a dispersive medium has been obtained from simple mathematical operations 
(derivatives) applied to the standard "X-wave". 



5. — Optical applications 

To illustrate what was said before, two practical examples will be considered, both 
in optical frequencies. When mentioning optics, it is natural to refer ourselves to optical 
fibers. Then, let us suppose the bulk of the dispersive medium under consideration to be 
fused Sihca (Si02). 

Far from the medium resonances (which is our case), the refractive index can be 
approximated by the well-known Sellmeier equation [18] 

nV)=l + E-#^, (22) 
j=i - ^ 

where uj is the resonance frequency, Bj is the strength of the jth resonance, and N is the 
total number of the material resonances that appear in the frequency range of interest. For 
typical frequencies of "long-haul transmission" in optics, it is necessary to choose N = 3, 
which leads us to the values[18] Bi = 0.6961663, B2 = 0.4079426, ^3 = 0.8974794, 
Ai = 0.0684043 fim, A2 = 0.1162414 fim and A3 = 9.896161 fim. 

Figure 7 illustrates the relation between N and uj, and specifies the range that will be 
adopted here. In the two examples, the spectra are localized around the angular frequency 
ujc = 23.56x10^^ Hz (which corresponds to the wavelength Ac = 0.8/im), with two different 
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bandwidth AcJi = O.SScJc and Aa;2 = OAuJc- The values of a and corresponding to 
these two situations are a = 1.14592 x 10~l^ = 27, and a = 1.90986 x 10"^^, = 41, 
respectively. 

Looking at these "windows", one can notice that Silica does not suffer strong vari- 
ations of its refractive index. As a matter of fact, a linear approximation to n = n{uj) 
is quite satisfactory in these cases. Moreover, for both situations, and for their respec- 
tive no values, the value of parameter 6 results to be very small, verifying condition (p!5[) : 
which means that it is quite acceptable our truncation of the Taylor expansion. The beam 
intensity profiles for both bandwidths are shown in Figs. 8 and 9. 

In the first figure, one can see a pattern similar to that of Fig. 2; but here, of course, 
the pulse is much more localized spatially and temporally (typically, it is a fentosecond 
pulse) . 

In the second figure, one can observe some little differences with respect to the first 
one, mainly in the spatial oscillations inside the wave envelope[19]. This may be explained 
by taking into account that, for certain values of the bandwidth, the carrier wavelength 
become shorter than the width of the spatial envelope; so that one meets a well defined 
carrier frequency. 

Let us point out that both these waves are transversally and longitudinally localized, 
and that, since the dependence of \I' on 2; and t is given hj ( = z — Vt, they are free from 
dispersion, just like a classical X-shaped wave. 

6. — Conclusions 

In this paper we have first worked out analytical Superluminal localized solutions to 
the wave equation for arbitrary frequencies and with adjustable bandwidth in vacum. The 
same methodology has been then used to obtain new, analytical expressions representing 
X-shaped waves (with arbitrary frequencies and adjustable bandwidth) which propagate 
in dispersive media. Such expressions have been obtained, on one hand, by adopting the 
appropriate spectrum (which made possible to us both choosing the carrier frequency 
rather freely, and controlling the spectral bandwidth), and, on the other hand, by hav- 
ing recourse to simple mathematics. Finally, we have illustrated some examples of our 
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approach with apphcations in optics, considering fused Sihca as the dispersive medium. 
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Figure Captions 



Figure 1 - Geometrical representation, in the plane {u!,kz), of the condition (2): see the 
text. 

Figure 2 - Illustration of the real part of an X-wave with bandwidth, Ao;, of 10 MHz 
starting from zero. 

Figure 3 - Behavior of the derivative number, m, as a function of the normalized band- 
width frequency, Au!±/u!c- Given a central frequency, cUc, and a bandwidth, Au!±, one 
finds the exact value of m by substituting these values into eq.(8.1). 

Figure 4 - Normahzed spectra for lOc — 23.56 x 10^^ Hz and different bandwidths. The 
first with N — 27 (solid line), and the second spectrum with N — Al (dotted hne). See 
the text. 

Figure 5 - The real part of an X-shaped beam for microwave frequencies in a dispersion- 
less medium. 

Figure 6 - Sketch of a generic device (axicon, hologram, etc.) suited to properly deviating 
the wave vector of each spectral component. 

Figure 7 - Variation of the refractive index n{u;) with frequency for fused Sihca. The solid 
hne is its behaviour, according to Sellmeir's formulae. The open circles and squares are 
the linear approximations ior N — 41 and N — 27, respectively. 

Figure 8 - The real part of an X-shapcd beam for optical frequencies in a dispersive 
medium, with N = 27. It refers to the larger window in Fig. 7. 

Figure 9 - The real part of an X-shaped beam for optical frequencies in a dispersive 
medium, with N = 41. It refers to the inner window in Fig. 7. 
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